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The  Dynamics  of  a  Rigid  Body  in  the  Space  Plasma 


1.  INTRODUCTION 

During  the  decade  of  the  1980s,  many  new  types  of  systems  will  be  placed  in 
orbit  around  the  earth.  They  will  have  a  "number  of  purposes,  the  two  most  prob¬ 
able  at  the  beginning  being:  facilities,  possibly  inhabited,  for  materials  produc¬ 
tion  and  processing  under  zero-gravity  conditions,  and  some  form  of  solar-power 
collector  with  a  means  of  transmitting  that  power  to  earth.  From  the  point  of 
view  of  the  design  and  on-orbit  operation  of  these  "Space  Stations,  "  the  most  dra¬ 
matic  change  from  the  orbiting  objects  of  the  1970s  will  be  their  size.  Early  ex¬ 
amples  of  these  systems  will  have  linear  dimensions  on  the  order  of  hundreds  of 
meters  up  to  several  kilometers,  with  sizes  projected  to  reach  10  km  by  the  mid- 
1990s.  1  With  the  exception  of  thin  antennas  and  booms  on  scientific  satellites, 
no  object  so  far  orbited  has  exceeded  a  few  tens  of  meters  maximum  dimension. 
The  structures  will  be  built,  initially,  by  deployment  of  assemblies  from  the 
Space  Shuttle  cargo  bay  and,  later,  by  making  use  of  the  multiple -flight  capability 
of  the  Space  Shuttle  to  transport  construction  materials,  fabrication  equipment, 
and  workers  to  a  low  earth  orbit. 

In  this  paper,  the  interaction  between  large  space  structures  and  their  j»nvi- 
ronment  will  be  considered  insofar  as  it  produces  mechanical  perturbing  forces 
(Received  for  publication  23  August  1979) 

1.  Hagler,  T. ,  Patterson,  H.  G. ,  and  Nathan,  C.  A.  (1977)  Learning  to  build 
large  structures  in  space.  Astronaut.  Aeronaut.  .  15:51. 
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that  might  impact  orbital  lifetime  or,  possibly  more  important  in  practical  terms, 
the  attitude  stability  and  control  of  such  a  structure.  We  will  not  consider  pos¬ 
sible  techniques  of  attitude  control,  but  simply  demonstrate  the  kinds  of  forces 
that  will  be  present.  The  structure  will  be  assumed  to  be  rigid,  non-rotating 
(see  Appendix  A),  and  any  oscillation  or  deformation  caused  by  the  external  forces 
will  not  be  considered  at  this  time.  It  should  be  remembered,  however,  that  the 
accuracy  and  reliability  of  attitude  control  is  crucial  in  some  applications,  such 
as  solar  power  transmission  to  earth.  In  those  cases  the  structures  will  have  to 
exhibit  rigidity  and/ or  controllability  over  a  wide  frequency  range  of  applied 

forces,  and  that  requirement  may  well  be  the  deciding  factor  in  the  choice  of  con- 

2 

struction  and  assembly  methods. 

There  are  many  independently  acting  physical  processes  that  act  upon  a  body, 
natural  or  artificial,  passing  through  the  space  plasma.  Frequently,  however, 
many  of  these  processes  result  in  mechanical  or  electromechanical  forces  on  the 
body  which  are  many  orders  of  magnitude  smaller  than  one  or  two  dominant  mech¬ 
anisms,  and  they  may,  therefore,  be  neglected.  Exactly  which  mechanism  will 
exert  the  greatest  force  on  a  body  in  a  specific  situation  depends  primarily  on  the 
body's  location,  size,  and  geometry.  Factors  such  as  solar  illumination  and 
spacecraft  attitude  may  cause  regular  o.  intermittent  changes  in  the  balance  of 
forces  acting  and  these  changes  can  have  a  cumulative  effect,  due  to  resonance 
with  the  structure  of  the  body  or  the  dynamics  of  its  orbit. 


2.  CLASSES  OF  INTERACTION  BETWEEN  A  MOVING  BODY  AND  A  PARTIALLY  IONIZED  PLASMA 

In  this  section,  the  "body"  (spacecraft  or  structure)  is  taken  to  be  electrically 
conducting,  non-magnetic,  and  to  have  no  rotation  about  its  center  of  mass  (see 
Appendix  A),  Unless  otherwise  stated,  no  specific  geometry  is  implied.  Colli¬ 
sions  between  particles,  charged  or  neutral,  are  neglected,  because  of  the  very 
low  collision  frequencies  above  about  150  km  altitude. 

2.1  Charging  Effects 

2.  1.  1  SPACECRAFT  POTENTIAL  (VEHICLE 
POTENTIAL  OR  PLASMA  POTENTIAL) 

Except  for  transients  caused  by  fluxes  of  energetic  charged  particles  or  by 
changing  solar  illumination,  there  can  be  no  net  flux  of  charge  to  an  isolated  body 
in  a  partially  ionized  plasma.  Different  parts  of  a  spacecraft  that  are  electrically 


2.  Oglevie,  R.  E.  (1978)  Attitude  control  of  large  solar  power  satellites,  AIAA 
Paper  78-1266,  p.  37,  in  Proceedings  of  AIAA  Guidance  and  Control  Con¬ 
ference.  Palo  Alto,  CA,  1978,  AlAA.  New  York. 
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isolated  can  charge  to  potentials  differing  by  thousands  of  volts.  This  is  the  case 
when  the  spacecraft  is  in  the  shadow  of  ihe  earth  at  geosynchronous  altitudes,  due 
to  the  low  ambient  thermal  plasma  density.  Electrical  discharge  between  surfaces 
can  occur  when  the  breakdown  voltage  is  exceeded.  If  this  happens  between  thin 
layers  of  dielectric  coatings  separating  electrical  components,  for  example,  it  can 
damage  or  destroy  critical  parts  of  electrical  systems.  It  is  the  need  to  eliminate 
this  damage  that  has  been  the  driving  force  behind  investigations  into  spacecraft 
charging  mechanisms  over  the  past  decade. 

The  effect  of  small  potentials  of  1  V  or  less  also  is  of  interest,  because  of  the 
drag  and  torque  forces  they  cause.  These  forces  are  caused  by  the  collision  and 
near  collision  of  the  charged  spacecraft  with  the  ambient  plasma  particles,  and 
also  by  the  interaction  between  the  ambient  magnetic  field  and  currents  flowing  to 
and  within  the  spacecraft. 

The  equilibrium  potential  $  of  a  body  in  the  space  plasma  is  the  potential  with 
respect  to  the  plasma  at  which  there  is  no  net  flow  of  charge  to  or  from  the  body. 
At  this  potential  <$,  we  have: 


=  l.  +  I  +  I  . 
+  S  ph 


(1) 


where  I_  is  incident  electron  current 
I+  is  incident  ion  current 

*s  is  a  positive  current  due  to  secondary  electrons  from  ion 
and  electron  impact  and  from  backscattering  of  incoming 
electrons 

Iph  is  photoemission  current. 

In  order  to  determine  4>,  expressions  must  be  developed  for  the  individual 
currents  constituting  Eq.  (1),  This  requires  data  on  the  thermal  energy  plasma, 
the  energetic  particle  fluxes,  and  on  the  backscattering  and  secondary  emission 
coefficients.  In  many  cases  these  data  are  not  available,  and  the  usual  practice 
is  to  use  laboratory  data  or  to  extrapolate  existing  space  measurements.  The  full 
version  of  Eq.  (1)  is  then  solved  numerically  for  $.  This  modeling  approach  has 
been  used  most  extensively  to  establish  the  equilibrium  potential  of  cylindrical-  or 
spherical -shaped  satellites  with  dimensions  of  a  few  meters.  In  such  models,  the 
relatively  small  disturbing  effect  of  the  earth's  magnetic  field  is  often  neglected; 
however,  the  fact  that  the  electrons  generally  can  travel  only  parallel  to  the 
magnetic  field  places  severe  limitations  on  possible  electron  fluxes  in  certain 

3.  Knott,  K.  (1972)  The  equilibrium  potential  of  a  magnetospheric  satellite  in  an 
eclipse  situation.  Planet.  Space  Sci.  .  20:1137. 


9 


geometries,  and  must  be  considered  in  any  detailed  treatment  of  the  equilibrium 
4  5 

potential.  ’  '  Magnetic  field  effects  will  be  considered  further  in  Sections  2.  5 
and  2.  G. 2. 

Under  most  conditions,  the  polarity  and  magnitude  of  $  depend  upon  the  bal¬ 
ance  in  Eq.  (1)  between  I  the  total  photocurrent,  and  the  portion  of  I_  that  is 
due  to  thermal -energy  electrons  reaching  the  body.  There  are,  however,  two 
important  exceptions  to  this  statement:  one  concerns  spacecraft  in  geosynchro¬ 
nous  orbit  under  eclipse  conditions  (spacecraft  in  the  earth's  shadow),  and  the 

other  concerns  regions  where  the  thermal  plasma  density  exceeds  a  value  of 
3  -3 

approximately  5  x  10  cm  .  In  the  first  case,  at  geosynchronous  altitude,  par¬ 
ticles  with  energies  on  the  order  of  several  tens  of  keV  frequently  are  present  and 
impinge  on  the  spacecraft.  If  electrons  predominate  in  this  flux,  then  in  shadow, 
where  no  photocurrent  is  generated  to  balance  this  negative  flux,  the  spacecraft 
will  rest  at  a  negative  potential  of  many  kilovolts  in  order  both  to  reduce  the  elec¬ 
tron  flux  and  to  attract  positive  ions  as  a  means  of  restoring  zero  net  charge  flow 

to  the  spacecraft.  In  the  second  case,  at  low  altitudes  in  a  relatively  dense  ther- 
3-3 

mal  plasma  (N  a  5  x  10  cm  ),  the  isotropic  flux  of  rapidly  moving  electrons  is 

-9 

several  orders  of  magnitude  greater  than  the  photocurrent  (i  ,  «  3  x  10  amp 
-  2  P** 

cm  ).  It  is  balanced  by  the  flux  of  relatively  slow  moving  positive  ions  swept  up 

by  the  rapid  motion  of  the  spacecraft  (vg  a  7  x  104 5 6 7 8  cm  sec  1),  and  under  these  con¬ 
ditions  the  absence  of  solar  illumination  to  generate  a  photocurrent  irmkes  little 
or  no  difference  to  the  value  of  the  equilibrium  potential  This  is  generally  the 
case  at  altitudes  below  about  1500  km  where  <S>  has  a  value  of  1  or  2  V  negative, 
although  plasma  densities  low  enough  to  allow  the  photocurrent  to  predominate  in 
Eq.  (1)  may  be  encountered  below  1500  km  in  auroral  and  polar  regions.  At 

_3 

plasma  densities  down  to  a  few  cm  ,  a  positive  potential  of  some  tens  of  volts  is 

sufficient  to  equalize  the  photocurrent  in  the  absence  of  high  energy  particles. 

As  indicated,  extensive  studies  have  been  made  of  the  equilibrium  conditions 

6  “  8 

of  small  spherical  satellites.  An  important  factor  in  those  studies  is  the  extent 


4.  Chu,  C.  K.,  and  Gross,  R.  A.  (1966)  Alfv6n  waves  and  induction  drag  on  long 

cylindrical  satellites,  AIAA  Journal.  4  (No.  12):2209. 

5.  Beard,  D.  B. ,  and  Johnson,  F.  S.  (1960)  Charge  and  magnetic  field  inter¬ 

action  with  satellites,  J.  Geophys.  Res.,  65:1. 

6.  Samir,  U.,  and  Willmore,  A.  P.  (1966)  The  equilibrium  potential  of  space¬ 

craft  in  the  ionosphere,  Planet.  Space  Sci.  ,  14:1131. 

7.  Chang,  H.  H.  C.,  and  Smith,  M.  C.  (1960)  On  the  drag  of  a  spherical  satel¬ 

lite  moving  in  a  partially  ionized  atmosphere,  Brit.  Interplanet.  Soc.  J.  , 
17:199. 

8.  Brundin,  C.  L.  (1963)  Effects  of  charged  particles  on  the  motion  of  an  earth 

satellite,  AIAA  Journal,  1  (No.  11):2529. 
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to  which  the  plasma  is  modified  by  the  electrostatic  charge  and  motion  of  the 
spacecraft.  This  falls  under  the  general  category  of  "sheaths  and  wakes"  phe¬ 
nomena,  and  is  of  great  importance  when  scientific  measurements  of  the  ambient 

pla.  na  are  required,  since  in  order  to  measure  charged  particles  with  energies 

9 

of  a  few  tenths  eV  the  collecting  probes  must  be  at  the  plasma  potential.  The 
electrostatic  charge  on  the  spacecraft  is  screened  from  the  ambient  plasma  by 
a  sheath  region  where  particles  of  one  sign  predominate  (positive  ions  in  the  case 
of  a  negatively  charged  spacecraft).  This  sheath  is  characterized  by  the  Debye 
length  (Xq),  which  is  the  distance  from  the  surface  at  which  the  potential  falls 
to  lie. 


(2a) 


(2b) 


where  Tg  is  thermal  electron  temperature  in  °K 

-3 

is  electron  density  in  n  cm 
e  is  electron  charge  in  Coulomb 

”16  ”1 

k  is  Boltzmann  constant  =  1.38  x  10  erg  °K 

The  degree  of  electrostatic  screening  and  the  effect  of  ion  and  electron 
thermal  velocities  in  "filling  in"  any  wake  region  left  by  the  rapid  passage  of 
the  spacecraft  through  the  plasma  must  be  considered  in  an  exact  evaluation  of 

o  in 

equilibrium  conditions.  '  Since  they  are  essentially  edge  effects,  however, 
sheaths  and  wakes  have  little  impact  on  the  mechanical  forces  acting  on  an  ex¬ 
tensive  planar  structure  in  the  space  plasma,  except  when  very  high  voltages  — 


9.  Samir,  U. ,  Maier,  E.  J. ,  and  Troy,  B,  E.  Jr.  (1973)  The  angular  distribu¬ 
tion  of  ion  flux  around  an  ionospheric  satellite,  J.  Atmos.  Terr. 

Phys.,  33: 513. 

10.  Kasha,  M.  A.  (1969)  The  Ionosphere  and  its  Interaction  With  Satellites. 

Chap.  3,  New  York,  Gordon  and  Breach. 
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on  the  order  of  thousands  of  volts— are  generated  (as  on  a  solar  power  station).  ’ 
We  note  here,  for  completeness,  the  equation  for  4  for  a  conducting  sphere, 
radius  R  in  the  150  km  to  1500  km  altitude  range.  At  such  altitudes  4  is  gener¬ 
ally  negative,  and  the  local  Debye  length  is  small  compared  with  the  dimensions 
of  the  sphere.  Neglecting  wake  effects  and  the  presence  of  the  earth's  magnetic 

O 

field,  4  is  given  by  Brundin  as: 


4 


-2  In 


1  - 


-£h 


v  tt  R  n  e  , 
s  e  > 


where  vfi  is  the  satellite  volocity 

m  is  electron  mass 
e 

I  ^  is  the  total  photoelectron  current. 


(3) 


2.  1.  2  PLANAR  SURFACE  CURRENT  FLOW 


When  considering  large-scale  conducting  structures,  it  is  necessary  to  con¬ 
sider  the  different  potentials  that  are  induced  at  different  locations  on  the  struc¬ 
ture  by  its  motion  through  the  earth's  magnetic  field  (Section  2.  5).  Currents  of 
different  polarities  and  area  densities,  therefore,  flow  through  the  surface  at 
different  locations  on  the  structure.  As  these  currents  flow  within  the  structure, 
they  cause  drag,  and  possibly  torque,  as  they  interact  with  the  ambient  magnetic 
field.  We  quote  here  expressions  for  the  currents  flowing  from  a  Maxwellian 

plasma  to  a  planar  surface  at  known  potential  for  both  accelerating  and  retarding 

1 3 

modes,  neglecting  edge  effects.  These  currents  must  be  evaluated  so  that  the 
flow  pattern  within  the  structure  can  be  used  to  calculate  the  resulting  magneti¬ 
cally  induced  drags  and  torques. 

11.  Parker,  L.  W.  (1979)  Plasma  sheath  effects  and  voltage  distributions  of 

large,  high -power  satellite  solar  arrays.  Proceedings  of  the  Spacecraft 
Charging  Technology  Conference,  USAF  Academy,  Colorado  Springs. 

Colo.  .  1978.  p.  341,  NASA  Conference  Publication  207 1 /AFGL-TR -79- 
0082. 

12.  Freeman,  J.  W.  ,  Cooke,  D.  ,  and  Reiff,  P.  (1979)  Space  environmental 

effects  and  the  solar  power  satellite.  Proceedings  of  the  Spacecraft  Charg¬ 
ing  Technology  Conference,  USAF  Academy,  Colorado  Springs,  Colo. , 
1978.  p.  408,  NASA  Conference  Publication  2071/ AFGL-TR -79-0082. 

13.  Whipple,  E.  C.  Jr.  (1959)  The  ion-trap  results  in  the  'exploration  of  the 

upper  atmosphere  with  the  help  of  the  third  Soviet  Sputnik,  '  Proc.  IRE. 
47:2023. 
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In  an  accelerating  mode  (♦  positive  for  electrons  and  negative  for  ions),  the 
current  flowing  to  area  A  is  given  by: 

i  -  rieA  ^ ,/n_ y  Cos  9  [  1  +  erf(y  Cos  8)  ]  +  1  exp  [  -  v2  Cos2  0  ]  I  . 

acc'v^)2  2  ('  “ 

In  a  retarding  mode  ($  negative  for  electrons  and  positive  for  ions). 


I  ,  =  neA 
ret  - 

\jv.  y  Cos  9  [  1  -  erf  (x-y  Cos  0  )  ]  +  1  exp  [  -(x-  y  Cos  8)2  ](  (5) 

1  9  9  ( 

\  tT3L 

I  2  J 

X  = 

v/e$7FF 

(6) 

a  = 

m/2kT 

(7) 

v  = 

(8) 

where  e  is  electron  charge 

n  is  density  of  the  particles 
m  is  mass  of  a  particle 

T  is  Maxwellian  temperature  of  the  particles 
k  is  Boltzmann  constant 

*  is  potential  difference  of  surface  relative  to  plasma 
A  is  area  of  the  surface 

vg  is  velocity  of  the  surface  relative  to  plasma 

0  is  angle  between  the  surface  normal  and  v  . 
b  s 


2.2  Solar  Radiation  Pressure 

Any  object  in  sunlight  will  be  subjected  to  a  furce  caused  by  momentum  trans¬ 
fer  from  solar  photons  to  any  surface  that  they  strike.  For  an  earth  satellite  in 
a  circular  orbit  that  does  not  enter  the  earth's  shadow,  the  net  accelerati  in  due 
to  solar  radiation  pressure  is  zero,  since  the  force  is  accelerating  for  cue  -half 
the  orbit  and  retarding  for  the  other  half.  For  a  body  in  an  elliptic  orbit  that 
enters  the  earth's  shadow,  there  will  be  a  net  change  in  the  energy  of  the  body 
due  to  solar  radiation  pressure,  unless  the  accelerating  and  retarding  radiation 
pressure  forces  are  equal  over  a  complete  orbit.  This  is  a  geometry  that  cannot 
persist  for  an  extended  period,  because  of  the  shift  of  the  orbit  in  inertial  space; 
therefore,  radiation  pressure  must  be  considered  in  the  calculation  of  most  or¬ 
bital  dynamics.  For  interplanetary  missions,  the  radiation  pressure  may  be  in 
the  same  sense  for  extended  periods  ar.d  can  be  harnessed  to  accelerate,  slow 


13 


14 

down,  or  stabilize  probes  during  long,  otherwise  unpowered,  coasting  phases. 

It  can  even  be  used  as  the  primary  means  of  propulsion,  as  in  the  NASA  solar 
sail  proposal  to  carry  a  payload  to  rendezvous  with  Halley's  comet. 

If  the  area  presented  to  solar  radiation  is  asymmetric  about  the  body's  center 
of  mass,  a  torque  will  be  applied  even  in  orbits  where  radiation  pressure  causes 
no  net  acceleration  of  the  center  of  mass.  This  mechanism  can  be  the  principal  cause 
of  applied  torque,  especially  in  the  case  of  a  structure  that  projects  large  optically 
opaque  surfaces  to  solar  radiation  (as  opposed  to  one  with  an  open  structure).  It 
would  be  possible,  under  some  orbital  and  spacecraft  configurations,  to  balance 
out  magnetic  and / or  aerodynamic  torques  by  careful  positioning  and  trimming  of 
open  and  closed  areas.  The  feasibility  of  using  solar  radiation  pressure  to  trim 
spacecraft  attitude  has  already  been  demonstrated  on  the  OTS-2  geosynchronous 
satellite.  (A  small  additional  force  under  conditions  of  solar  illumination 
arises  from  the  recoil  of  the  solar-induced  photoelectrons  from  the  body.  This 
is  on  the  order  of  1  percent  of  the  true  radiation  pressure  and  may  be  neglected. ) 

For  a  perfectly  absorbing  surface,  the  radiation  pressure  is  equal  to  the 
energy  density  of  solar  radiation,  and  the  force  acting  on  a  body  is  given  by: 


F„  =  L  A  Cos  a  dyne,  (9a) 

K  _ 

4rrRs  c 

33  -2  -I 

where  L  is  total  energy  radiated  from  the  sun  (3.86  x  10  erg  cm  sec  ) 

R  is  distance  from  the  sun  in  cm 

s  10-1 
c  is  velocity  of  light  (3  x  10  cm  sec  ) 

A  is  area  exposed  to  solar  radiation 

a  is  angle  between  the  normal  to  A  and  the  sun-body  line. 

In  the  neighborhood  of  the  earth,  Rg  =  1  AU  (1.495  x  101;i  cm)  and  the  radia¬ 
tion  pressure  Cr  is  given  by: 


C  r  4.6  x  10  Cos  or  dyne  cm 


-2 


(9b) 


14.  Sohn,  M.  (1959)  Attitude  stabilization  by  means  of  solar  radiation  pressure. 

Am.  Rocket  Soc.  J.  29:371. 

15.  Renner,  U.  (1979)  Attitude  control  by  solar  sailing  -  a  promising  experiment 

with  OTS-2,  European  Space  Agency  Journal  3:35. 
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2.3  Neutral  Particle  lira' 


Up  low  about  400  km  altitude,  the  largest  contribution  to  spacecraft  decelera¬ 
tion  comes  from  collisions  with  the  neutral  constituents  of  the  atmosphere,  which 
greatly  outnumber  the  charged  particles  in  these  regions.  Conditions  at  all  prac¬ 
tical  altitudes  for  spacecraft,  above  about  150  km,  are  those  of  "free  molecular 
flow,  "  meaning  that  the  mean  free  patli  of  all  particles  is  much  greater  titan  the 
dimensions  of  the  body  with  which  they  interact.  Additionally,  for  regions  where 
neutral  particle  collisions  are  significant,  their  thermal  velocities  arc  much 
lower  than  the  velocity  of  the  spacecraft  (vg  7x10'^  cm  sec  *)  .  The  magnitude 
of  the  drag  force  under  these  conditions  is  given  by: 


FND  t  I  P V 


<dA 


dvne, 


(10) 


where  l'  ^D  is  the  drag  force 

p  is  the  neutral  atmosphere  density  in  gm  cm  '* 
v  is  the  spacecraft  velocity  in  cm  sec 
Cn  is  the  drag  coefficient 

u  2 
A  is  cross-seetional  area  projected  into  the  flow  in  cm  . 

The  value  of  depends  upon  the  geometry  of  the  body,  the  mode  of  reflection  of 
the  incoming  particles  (specular  or  diffuse),  and  the  degree  of  thermal  accomoda¬ 
tion  of  the  particles  to  the  surface  temperature  of  the  body.  A  value  of  2.  2  is 
generally  used  for  for  spherical  and  nearly  spherical  bodies  in  the  150  to 
400  km  altitude  range.  Shapes  such  as  flat  plates  and  cylinders  have  somewhat 

1  fi 

higher  values  of  C^,  although  this  is  also  dependent  on  the  reflection  processes. 

At  higher  altitudes,  above  about  GOO  km,  the  thermal  velocity  of  all  particles, 

neutral  and  charged,  must  be  considered;  *  1  Cn  rises  to  a  value  of  greater  than 

^  16  18 

4.0,  and  the  drag  process  becomes  more  efficient.  ’  Because  of  the  much 
greater  degree  of  ionization  at  these  higher  altitudes,  however,  the  relative  con¬ 
tribution  of  neutral  particle  collisions  to  the  total  drag  becomes  very  small. 


16.  Cook,  G.  E,  (1965)  Satellite  drag  coefficients.  Planet.  Space  Sci.  13:929. 

17.  Schamberg,  R.  (1 959)  Analytic  representation  of  surface  interaction  for  free 

molecule  flow  with  application  to  drag  of  various  bodies,  in  Aerodynamics 
of  the  Upper  Atmosphere,  edited  by  D.  J.  Masson,  p.  12-1,  Rand 
Corporation. 

18.  Epstein,  P.  S.  (1924)  Resistance  experienced  by  spheres  in  their  motion 

through  gases,  Phvs.  Rev.  23:710. 
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2.4  Charged  Particle  Drag  (Coulomb  Drag) 


2.4.  1  DIRECT  IMPACT  DRAG 

The  potential  of  the  spacecraft  modifies  the  trajectories  of  charged  particles 
that  otherwise  would  not  strike  it.  Since  electrons  may  be  neglected  in  drag  cal¬ 
culations,  due  to  their  low  mass  compared  with  any  ion,  the  usual  situation  is 
that  the  effective  collisional  area  of  the  spacecraft  is  increased  by  the  attraction 
of  positive  ions  to  a  negatively  charged  body,  and  decreased  by  a  positively 
charged  body. 

In  the  simplest  case,  for  a  spherical  body  moving  at  VgS&Vj  where  v^  is  the 
ion  thermal  velocity,  and  when  all  ions  stick  to  the  body  after  impact  and  neutrali¬ 
zation,  the  drag  is  simply: 


Fit,  =  1  ~  m.n.v 

ID  2  its 


dyne. 


(11) 


where  m.  is  ion  mass  in  gm 

1  -3 

n.  is  ion  density  in  n  cm 

1  _i 

vg  is  the  satellite  velocity  in  cm  sec 

r  is  effective  radius  for  ion  collection  in  cm. 
e 

From  consideration  of  momentum  conservation,  rg  is  given  by: 

1 


re=R 


1  +  2e$ 


(12) 


where  R  is  the  radius  of  the  spherical  body  and  #  is  the  spacecraft  potential. 

The  exact  determination  of  the  effect  of  spacecraft  charge  depends  upon  the 
degree  to  which  the  charge  is  screened  from  the  surrounding  plasma.  This  is 
usually  defined  in  terms  of  the  local  Debye  length  (lip),  but  it  must  be  remembered 
that  this  is  the  distance  at  which  the  potential  due  to  spacecraft  charge  is  reduced 
by  l/e,  and  it  is  not  the  limit  of  the  disturbed  region.  Where  Xp  is  more  than  a  few 
centimeters,  the  disturbed  region  around  the  spacecraft  can  be  on  the  order  of 
10  to  40  percent  of  the  total  cross  section  of  a  spacecraft  1  or  2  m  in  diameter.  ® 

At  impact,  the  ions  are  neutralized  and,  if  re-emitted,  carry  momentum  that 

17 

is  a  further  source  of  drag  if  the  reflection  mechanism  is  diffuse.  This  re¬ 
emission  drag  can  be  greater  than  the  initial  impact  drag  of  the  ions  if  there  is  no 


16 


thermal  accommodation  of  the  incident  particles  to  the  body  surface  temperature. 

O 

This  drag  has  a  maximum  value,  F  ,  given  by  Brundin  as: 


I  dyne. 

i 

where  n.  is  ion  number  density  in  n  cm'J 
m.  is  ion  mass  in  gm 
R  is  the  spherical  body  radius  in  cm 
♦  is  spacecraft  potential  in  volts 
vg  is  the  satellite  velocity  in  cm  sec  *  . 


(13) 


2.4.2  NON-COLLISIONAL  DRAG  (MAXWELL 
DRAG  OR  DYNAMIC  FRICTION) 

An  additional  source  of  charged  particle  drag  arises  from  Coulomb  interaction 
leading  to  momentum  exchange  between  a  charged  body  and  ions  that  do  not  strike 
it  but  whose  trajectories  are  changed.  A  critical  factor  in  determining  this  drag 
is,  again,  the  extent  to  which  the  effect  of  spacecraft  potential  penetrates  into  the 
ambient  plasma.  The  problem  has  been  approached  by  a  number  of  workers7, 
who  made  varying  assumptions  about  the  extent  and  geometry  of  the  electrostatic 
screening  of  a  body  moving  through  a  plasma;  however,  in  considering  large  struc¬ 
tures,  we  need  note  only  that  non-collisional  drag  is  a  significant  fraction  of  the 
total  charged -particle  drag  only  when  the  body  in  question  has  dimensions  no  more 
than  an  order  of  magnitude  greater  than  the  local  Debye  length.  In  regions  where 
this  condition  is  satisfied,  for  very  large  structures  in  geosynchronous  orbit  out- 

_3 

side  the  plasmapause,  for  example,  ion  densities  are  too  low  (c  10  cm  )  to  con¬ 
tribute  any  significant  drag. 


2.3  Magnetic  Field  Drag  Effects 

The  effects  of  the  earth's  magnetic  field  on  the  potential  distribution  and  cur¬ 
rent  flows  on  a  conducting  body  in  the  space  plasma  are  not  significant  in  practice 

19.  Jastrow,  R. ,  and  Pearse,  C.  A.  (1957)  Atmospheric  drag  on  the  satellite, 

J.  Geophys.  Res.  62:413. 

20.  Fournier,  G.  (1970)  Electric  drag,  Planet.  Space  Sci.  18:1035. 

21.  Knechtel,  E.  D.  ,  and  Pitts,  W.  C.  (1964)  Experimental  investigation  of 

electric  drag  on  satellites,  AIA A  Journal,  2  (No.  6):1148. 
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unless  the  dimensions  of  the  body  exceed  approximately  10  m.  As  linear  dimen¬ 
sions  increase,  the  potentials  generated  on  the  body  by  its  motion  across  the  mag¬ 
netic  field  begin  to  exceed  the  energies  of  the  plasma  particles  with  respect  to  the 
body.  These  usually  lie  within  the  range  of  ±  25  eV.  Ions  and  electrons  flow  pref¬ 
erentially  to  the  negative  and  positive  sections  of  the  body;  therefore,  although 
there  is  a  zero  net  current  flow  to  the  whole  body,  positive  and  negative  currents 
flow  to  and  within  specific  parts  of  the  body.  This  results  in  drag  (Lorentz  forces) 
on  the  body  as  the  various  currents  interact  with  the  ambient  magnetic  field. 

For  a  conducting  body,  as  shown  in  Figure  1,  with  length  2L^  perpendicular 
to  vg  and  B,  a  potential  Ex  will  be  set  up  along  the  length  of  the  body: 

=  10  ^(vg  x  B)  volts  cm  * ,  (14) 


v  is  in  cm  sec 
s 

B  is  in  gauss. 

Figure  2  shows  the  distribution  of  the  currents  that  will  flow,  partly  as  a  result  of 
this  potential  and  partly  due  to  collisions  with  ambient  ions  and  electrons,  when 
L  is  on  the  order  of  a  few  meters.  The  induction  drag  or  Lorentz  force  (F,  ) 
acting  on  the  body  as  a  result  is: 


fl  =  V 

x 


I(x)sin  6  dx  dyne. 


where  I(x)  is  the  net  current  flowing  in  an  increment  dx  of  the  body  located  at  x 
6  is  the  angle  between  the  direction  of  current  flow  and  the  magnetic  field  B. 
For  the  configuration  of  Figure  1,  6  =  90 J  ,  and  if  the  currents  are  equally 
distributed  about  the  center  of  the  structure,  they  flow  for  an  average  distance  of 
L  perpendicular  to  B,  lienee: 


i'l  =10  IB  1,  dyne. 


where  1  is  the  mean  current  in  amps 
B  is  in  gauss 


1.  is  in  cm. 


Figure  1.  The  Basie  Structure  Moving  Perpen¬ 
dicular  to  the  Local  Magnetic  Field 


Figure  2.  The  Current  and  Potential  Distri 
butions  Along  the  x  Axis  of  the  Structure, 

When  [,  is  on  the  Order  of  a  Few  Meters 
x 

1!* 


The  location  of  "a,  "  the  point  on  the  structure  that  rests  at  plasma  potential, 
is  determined  by  balancing  the  electron  flux  to  the  positive  part  of  the  structure 
with  the  ion  flux  to  the  negative  part  and  the  photocurrent  (positive)  emitted  pre¬ 
dominantly  from  the  negative  section.  Since  the  electron  flow  to  the  structure  is 
parallel  to  the  magnetic  field,  current  densities  and  distributions  depend  on  the 
orientation  of  the  structure. 

3  -3 

At  altitudes  where  plasma  densities  are  approximately  5  x  10  cm  or 
greater,  point  "a,  "  for  structures  no  more  than  a  few  meters  in  size,  is  generally 
located  between  one-third  and  one-half  of  the  length  of  the  structure  from  the  posi¬ 
tive  end.  As  the  linear  size  is  increased  along  the  x  axis,  higher  potentials  are 

induced  across  the  structure  and  the  positive  end  is  effectively  grounded  in  the 

4 

local  plasma;  that  is  to  say,  a  »  L^.  At  higher  altitudes  for  structures  of  all 
sizes,  the  photocurrent  is  significantly  greater  than  either  the  ion  or  electron  flux. 
The  entire  strucutre  rests  at  a  positive  potential  with  respect  to  the  plasma  and 
attracts  an  increased  flux  of  electrons  to  balance  the  photocurrent.  Electrons  still 
flow  preferentially  to  the  less  negative  end  of  the  structure,  and  the  magnetically 
induced  potential  across  the  structure  will  still  be  generated. 


2.6  Torque  Processes 

All  the  drag  processes  outlined  so  far  also  will  exert  a  torque  on  the  body  if 
the  forces  produced  have  a  non-zero  total  moment  about  the  body's  center  of  mass 
Torque  also  will  be  generated  in  the  presence  of  a  magnetic  field  by  current  loops 
within  the  body,  but,  in  this  case,  no  additional  drag  results  as  long  as  the  mag¬ 
netic  field  is  uniform  over  the  dimensions  of  the  body.  This  process  may  have 
serious  consequences  when  dealing  with  large  structures  such  as  solar  power  sta¬ 
tions,  where  very  large  currents  are  generated.  Careful  routing  of  current  paths 
is  required  to  minimize  the  total  magnetic  moment  of  the  structure. 


2.6.  1  MOMENTUM  TRANSFER  TORQUES  (PARTICLE 

IMPACTS  AND  RADIATION  PRESSURE  TORQUES) 

If  a  drag  process  results  in  a  force  P  per  unit  area  projected  in  the  direction 
of  the  body's  velocity  vector  vg,  then  the  total  torque  T^  is: 


Tp  -  Pj*j*  r  dA  Cos  9  dr  dyne  cm. 


(17) 


whore  P  is  drag  pressure  in  dynes  cm 

r  is  vector  distance  of  dA,  from  center  of  mass  in  cm 

2 

dA  is  element  of  surface  area  in  cm 

v  is  %'elocitv  in  cm  sec  * 
s 

8  is  the  angle  between  v  and  the  normal  to  dA. 


2.6.2  MAGNETIC  TORQUES  (INDUCTION  TORQUES) 

When  induction  drag  results  from  currents  flowing  through  the  body  from  the 
space  plasma,  a  torque  T^  will  result  if  there  is  a  net  non-zero  moment  of  these 
drag  forces  about  the  center  of  mass.  From  Eq.  (15)  we  obtain: 


TL  10  B'lx  J  X  I(x)  Lx  Sin  5  dx  dyne  cm. 


where  B  is  magnetic  field  in  gauss 

I(x)  is  current  in  amps  flowing  in  an  element  dx  located  at  Lx  cm  from  the 
center  of  mass 

6  is  the  angle  between  the  direction  of  current  flow  and  the  magnetic 
field. 

In  addition,  a  current  flowing  around  a  closed  loop  in  the  presence  of  a  uni¬ 
form  magnetic  field  generates  a  torque  given  by: 


T.  =  10~  1  I  A  B  Sin  a'  dyne  cm. 


where  I  is  current  around  loop  in  amps 

2 

A  is  area  enclosed  by  loop  in  cm 
B  is  magnetic  field  in  gauss 

a'  is  angle  between  the  magnetic  field  and  the  normal  to  the  plane  of  the 
loop. 


2.  6.  3  GRAVITATIONAL  TORQUE 


The  gravitational  potential  gradient  of  the  earth  exerts  a  torque  on  a  body  in 
such  a  way  that  its  axis  of  minimum  moment  of  inertia  aligns  itself  with  the  radius 
through  the  center  of  the  earth;  that  is  to  say,  the  long  axis  of  a  body  will  line  up 
with  the  local  vertical.  Use  is  made  of  this  torque  in  the  technique  of  "gravity 


gradient  stabilization''  for  satellites  that  are  required  to  present  the  same  side 

toward  the  earth  at  all  times.  Stabilization  is  achieved  by  extending  a  mass  at 

the  end  of  a  long  boom  erected  along  the  axis  required  to  point  downwards.  This 

greatly  increases  the  moments  of  inertia  about  the  axes  perpendicular  to  the  boom 

22 

which  then  aligns  itself  along  the  local  vertical. 

Gravitational  torques  have  had  little  effect  on  the  stability  and  alignment  of 
many  satellites  launched  to  date  because  of  their  relatively  small  overall  dimen¬ 
sions,  and  the  fact  that  their  moments  of  inertia  about  all  three  axes  are  fre¬ 
quently  of  the  same  order  of  magnitude.  In  addition,  antennas  and  solar-cell 
arrays  have  made  large  relative  contributions  to  the  cross-sectional  area  of  the 
spacecraft  compared  with  small  additions  to  moments  of  inertia.  Because  of  this, 
particle  impacts,  magnetic  induction  effects,  and,  at  higher  altitudes,  solar  radia 
tion  pressure,  have  been  the  principal  sources  of  torque.  This  situation  is 
changed,  however,  when  we  consider  massive,  elongated  structures  with  high 
transparencies  for  particle  collision  and  solar  radiation  pressure. 

Figure  3  shows  the  previous  sample  structure  aligned  at  an  angle  0  to  the 
local  vertical  r  such  that  x,  7,  and  r  are  coplanar.  The  gravitational  torque  T 
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about  the  y  axis  is  given  by  Thomson  as: 

T  "  3GM  [  (Ml)  -  (Ml)  ]  Sin  2  3  dyne  cm,  (20) 

jo  z  x 

2n 

-8  2  ~2 

where  G  is  universal  gravitation  constant  =  6.  67  x  10  dyne  cm  gm 
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M  is  mass  of  the  earth  =  5.  973  x  10  gm 

R  is  distance  to  center  of  the  earth  in  cm 

2 

(Ml)x  are  moments  of  inertia  about  x  and  z  axes  in  gm  cm 
3  is  the  angle  between  the  local  vertical  r,  and  the  x  axis  of  the 

structure. 


22.  Fischell,  R.  E.  ,  and  Mobley,  F.  F.  (1964)  A  system  for  passive  gravity  - 

gradient  stabilization  of  earth  satellites,  in  Progress  in  Astronautics  and 
Aeronautics:  Guidance  and  Control  -  II.  Vol.  13,  p.  37,  New  York, 
Academic  Press. 

23.  Thomson,  W.  T.  (1962)  Spin  stabilization  of  attitude  against  gravity  torque, 

J.  Astronaut.  Sci.  9:31. 
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GRAVITATIONAL  TORQUE: 

^[(mI)z-(mI)x]Sin2  0 

Figure  3.  Geometry  to  Produce  a  Gravi¬ 
tational  Torque  About  the  y  Axis  of  the 
Structure.  The  x  and  z  axes  are  coplanar 
with  r,  the  local  vertical 


3.  NUMERICAL  EXAMPLES  -  DRAG 

In  this  section  we  calculate  values  for  the  drag  forces  on  a  sample  structure 
that  are  generated  at  a  number  of  altitudes  by  the  processes  outlined  in  Section  2. 
Representative  values  of  environmental  parameters  are  given  in  Table  1  for  an 
exospheric  temperature  of  1500°K.  The  structure  is  as  shown  in  Figure  1  and  its 
dimensions,  mass,  and  velocity  at  the  various  altitudes  are  given  in  Table  2. 

Table  1.  Environmental  Parameters 


Altitude 
(h  km) 

250 

500 

1,000 

36,  000(6.  6Re) 

Neutral  Density 
(P  gm  cm"3) 

10'13 

3xl0‘15 

2xl0'13 

— 

Neutral  Density 
(N  cm*3) 

109 

108 

105 

— 

Plasma  Density 
(m  cm"3) 

5xl05 

106 

3xl04 

5x10°  -  5xl02 

Neutral  Temperature 
(Ti«Tn  °K) 

800 

1,000 

2,000 

5,000  -  20,000 

Electron  Temperature 
(Te°K) 

1,000 

1,  500 

2,000 

5,000  -  20,000+ 

(inside  plasmasphere) 
Several  keV 
(outside  plasmasphere) 

Mean  Ion  Mass  (m  amu) 

24 

16 

8 

1 

Debye  length  (X^cm) 

0.3 

1.0 

2.0 

200  + 

Geomagnetic  Field 
(B  gauss) 

0.45 

0.40 

0.  33 

10'3 

Table  2.  Structure  Details 


Dimensions  (cm):  Length 

(21.  ) 

2  x  10  :’ 

X 

Height 

(2Lz) 

103 

Thickness 

(2V 

.102 

Total  Projected  Area 

(2Lx  x  2Lz  cm2) 

2  x  108 

Mass  per  Unit  Projected  Area 

(gm  cm  2) 

101 

Total  Mass  (gm) 

0 

2  x  10J 

Velocity  (cm  sec  8): 

250  km 

7.  80  x  10f> 

500  km 

7.60  x  lu5 

1 , 000  km 

7. 35  x  10° 

36,  000  km 

3.00  x  10S 

3.1  Solar  Radiation  Pressure 

In  the  near-earth  region  (1  All),  the  value  of  solar  radiation  pressure  is  ef¬ 
fectively  constant  at  all  altitudes.  For  a  totally  absorbing  surface,  the  pressure 
is: 


C  =4.6x10  Cos  a  dyne  cm  2  (9a) 

If  the  structure  faces  directly  toward  the  sun,  a  -  0”  and  the  total  force  acting 
8  2 

on  an  area  of  2  x  10  cm  will  be: 

4.  6  x  10'ax2x  108  =  9.  2  x  103  dyne. 
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3.2  Neul ral  Particle  Drag 


The  neutral  particle  drag  is  a  decelerating  force,  acting  at  all  times  in  the 
opposite  direction  to  the  velocity  vector  vg: 


FND 


A. 


(10) 


From  Section  2,3,  the  Drag  Coefficient,  C„,  =  2.2,  and  the  area  projected  into  the 

8  2  u 

velocity  vector.  A,  =2x10  cm  .  The  neutral  drag  at  three  altitudes  is  shown  in 
Table  3.  Neutral  densities  are  considered  negligible  at  geosynchronous  orbit.  Note 
that  at  an  altitude  of  1000  km  the  neutral  component  is  already  less  than  the  solar 
radiation  pressure  drag  (Section  3.  1)  at  the  time  when  the  strucutre  is  facing  directly 
into  the  sun;  the  net  energy  transfer  due  to  solar  radiation  pressure,  however,  is 
very  close  to  zero  over  a  complete  orbit. 


Table  3.  Neutral  Drag 


h  (km) 

250 

500 

1000 

p  (gm  cm  3) 

10'13 

3  x  1015 

2  x  10'17 

v  (cm  sec  *) 
s 

7.8  x  105 

7.6  x  105 

7.  35  x  105 

FND  (dyne) 

1.  34  x  107 

3.81 x  105 

2.38  x  103 

_  2 

Fnd/A  (dyne  cm  ) 

6. 69  x  10~2 

1. 91  x  10"3 

1. 19  x  10~5 

3.3  Charged  Particle  Drag 


The  charged  particle  drag  is  a  decelerating  force  acting  in  the  same  sense  as 
neutral  particle  drag: 


FID  =lmiVs2A 


dyne. 


(21) 


where  rm  is  mean  mass  of  ambient  ions  in  gm 

n.  is  number  density  of  ambient  ions  in  n  cm 

v  is  the  spacecraft  velocity  in  cm  sec  ^ 

s  2 
A  is  the  area  projected  into  vg  in  cm  . 
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The  values  of  drag  given  in  Table  4  are  derived  from  Eq.  (21)  for  a  structure 

8  2 

of  cross-sectional  area  A  =  2  x  10  cm  moving  much  faster  than  the  thermal 
velocity  of  the  ions.  No  allowance  is  made  for  increased  effective  area  due  to 
Coulomb  attraction  of  nearby  ions,  since  that  is  most  significant  when  the  dimen¬ 
sions  of  the  structure  are  no  more  than  one  order  of  magnitude  greater  than  the 
local  Debye  length.  Similarly,  no  diffuse  re-omission  of  ttie  neutralized  ions  is 
considered.  Since  both  of  these  factors  would  increase  charged  particle  drag,  the 
figures  given  in  Table  4  are  conservative.  Since  thermal  ion  velocities  at  geo¬ 
synchronous  orbit  are  greater  than  orbital  velocity,  Eq.  (21)  does  not  apply;  how¬ 
ever,  the  ambient  density  is  so  low  at  this  altitude  that  the  drag  resulting  from 
ion  impact  is  negligible  compared  with  solar  radiation  pressure  forces. 


Table  4.  Charged  Particle  Drag 


h  (km) 

250 

500 

1000 

m.  (gm) 

4.01  x  10'23 

_9? 

2.68  x  10 

.  ,n-23 

1. 34  x  10 

n.  (n  cm  3) 

5  x  105 

106 

3  x  104 

vg  (cm  sec  1) 

7.8  x  105 

7. 6  x  105 

7. 35  x  105 

fid  (dyne) 

1.  22  x  103 

1.  54  x  103 

2.  17  x  101 

_  2 

Fjd/ A  (dyne  cm  ) 

6.  1  x  10'6 

7.73  x  10~6 

1.08  x  10^7 

3.4  Magnetic  Drag  (Induction  Drag) 

The  three  processes  dealt  with  so  far  all  depend  on  the  area  that  the  structure 
projects  into  the  velocity  vector  or  the  sun  vector.  If  an  open-girder  type  of  con¬ 
struction  is  used,  the  forces  resulting  from  these  processes  would  be  reduced  in 
direct  proportion  to  the  optical  transparency  of  the  area  projected.  This  change 
in  opti  :al  transparency  must  be  uniform  ever  the  whole  structure,  however,  or  a 
torque  will  result. 

I  i  the  case  of  magnetic  field  drag  (induced  drag),  it  is  the  linear  dimensions 
and  orientation  of  the  structure  in  the  local  magnetic  field  that  are  critical.  From 
Eq.  (16),  the  drag  for  the  orientation  in  Figure  1  is: 

F,  =  10'1  I  B  L  dyne.  (16) 

Li  X 
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In  practice,  the  actual  value  of  I  depends  upon  the  local  ion  density,  solar  illumi¬ 
nation,  transparency,  the  velocity  of  the  structure,  and  on  orientation  in  the  mag¬ 
netic  field  as  it  affects  the  potential  distribution  over  the  structure.  For  this  ex¬ 
ample  we  consider  only  the  photocurrent  (Ip)  due  to  solar  UV  photons  and  the  ion 
collection  current  (I.)  swept  up  by  the  structure  whose  velocity  is  assumed  to  be 
large  compared  with  the  ion  thermal  velocities.  This  assumption  is  not  true  at 
geosynchronous  orbit  altitudes  where  the  photocurrent  dominates.  In  an  equato¬ 
rial  orbit,  the  photocurrent  contribution  will  be  zero  at  1200  hr  I.T,  and  maximum 
at  0G00  hr  LT  and  1800  hr  LT,  if  both  back  and  front  surfaces  are  taken  to  be 
photoemitters.  If  almost  all  of  the  surface  of  the  structure  is  an  ion  collector 
and/or  a  photoemitter,  then  the  average  distance  perpendicular  to  B  through  which 
these  currents  flow  will  be  close  to  (  =  10  cm),  one-half  the  total  length.  At 
each  altitude  we  find  I.,  using  the  environmental  data  in  Table  1,  and  then  substi- 

1  -  j 

tute  in  Eq.  (16)  to  determine  F,  .  Photocurrent  I  -  6  x  10  amp  at  all  altitudes, 

h  P 

and  its  variation  around  the  orbit  is  reflected  in  the  range  of  values  for  I.  given 
in  Table  5. 


Table  5.  Induction  Drag 


h  (km) 

250 

500 

1000 

36000 

I.  (amp) 

(1.  25-1.  3DX10'1 

(2.43-2.49)xl01 

(T.oe-is.oeixio'1 

6x10  1 

B  (gauss) 

4. 5x10* 

4.0xl0_1 

3.3x10'* 

10"3 

F,  Max 

5.  9xl04 

9.96X104 

4. 32xl03 

6x10° 

F.  /A 

2.  94x10  ^ 

4.98xl0'4 

2. 16xl0'5 

3xl0"8 

(dyne  cm 

4.  NUMERICAL  EXAMPLES  -  TORQUE 

4.1  Momentum  Transfer  Torques 

In  order  to  estimate  the  magnitude  of  these  torques  (due  to  solar  radiation 
pressure  and  collisions  with  charged  and  neutral  particles),  we  modify  the  struc¬ 
ture  as  shown  in  Figure  4.  One-quarter  of  the  length  is  covered  with  a  conductor, 
as  before,  but  the  remainder  now  has  an  open  construction  with  an  optical  trans¬ 
parency  of  90  percent.  The  total  mass  and  its  distribution  about  the  center  of 
gravity  remains  unchanged.  The  y  axis  of  the  structure  is  maintained  parallel  to 

v  at  all  times, 
s 
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OPTICAL  TRANSPARENCY 

SECTION  1=  90% 
SECTION  0% 


Figure  4.  Modification  of  the  Basic  Structure  to- 
Demonstrate  Momentum  Transfer  and  Solar  Radia¬ 
tion  Pressure  Torques  by  the  Imbalance  of  Forces 
About  the  Center  of  Gravity 

The  applied  torque  will  have  a  maximum  value  when  the  structure  is  moving 
directly  towards  the  sun.  The  solar  radiation  pressure  component  varies  con¬ 
tinuously  throughout  the  orbit  and  acts  to  equalize  the  particle  impact  torques 

whenever  v  has  a  component  away  from  the  sun. 

s  -2 
If  the  pressure  on  the  structure  at  any  time  is  P  dyne  cm  ,  then  from 

Eq.  (17)  the  resulting  torque  Tp  will  be: 


+L. 


T=2LP  f  xL  dl.  d'ne  cm. 
p  z  T  xx- 


(22a) 


-L 


Due  to  the  basic  symmetry  about  the  z  axis,  Tp  arises  from  the  additional 
area  in  the  enclosed  section  of  the  structure  so  that: 


T  =0.9',  2L  P 
P  |  z 


L  2  L  2 

X  -  X 


J) 


(22b) 


T  =  3.38  x  10 *2  P  dyne  cm. 
P 


From  Section  3  we  obtain  maximum  values  for  P  in  Table  6  by  adding  solar 
radiation  pressure  and  the  neutral  and  charged  particle  forces. 


i 

i 
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Table  G.  Pressure  Torque 


h  (km) 

250 

500 

1000 

36000 

-2 

P(dyne  cm  ) 

6.69  x  10"2 

1.  91  x  10'3 

5.80  x  10'5 

4.  6  x  10'5 

Tp(dyne  cm) 

2.  26  x  1011 

6.46  x  109 

1.  96  x  108 

1.  56  x  108 

4.2  Magnetic  Torques  (Induction  Torques) 

As  an  example  of  currents  flowing  in  a  section  of  the  structure  to  produce  a 
torque  by  imbalance  about  the  center  of  mass,  we  will  use  the  currents  derived 
in  Section  3.4.  We  reduce  them  by  a  factor  of  4  to  represent  their  flowing  in  the 
enclosed  portion  of  the  structure,  and  ignore  currents  flowing  in  the  rest  of  the 
structure.  The  mean  moment  arm  of  the  Lorentz  forces  produced  about  the  cen¬ 
ter  of  mass  will  be  3L  /4,  and  the  currents  will  flow  for  a  mean  distance  of  L  /4 
x  x 

perpendicular  to  B.  In  Table  7  we  use  these  values  in  L'q.  (18)  to  obtain  T^. 


Table  7.  Induction  Torque 


h  (km)  250 

500 

1000 

36000 

I  (amp)  3.  28  x  10° 

6.  23  x  10° 

3.27  x  10'1 

1. 50  x  10'1 

B  (gauss)  4.  5  x  10* 1 

4.0  x  10*1 

3.  3  x  10_1 

IQ'3 

Tl  (dyne  cm)  2.  77  x  108 

4.67  x  108 

2.03  x  107 

2.81  x  104 

To  estimate  the  torque  (T.)  caused  by  a  flow  of  current  around  a  closed  loop 
within  the  structure,  we  postulate  a  net  current  of  1  amp  flowing  around  the  outer¬ 
most  edges  of  the  structure,  aligned  as  in  Figure  1,  so  that  the  torque  acts  about 
the  x  axis  of  the  structure.  Note  that  for  this  source  of  torque  internal  construc¬ 
tion  and  optical  transparency  are  irrelevant;  the  torque  is  determined  by  the  area 
enclosed  by  the  loop  and  its  alignment  in  the  magnetic  field.  In  Table  8  we  use 
Eq.  (19)  and  Table  1  to  obtain  values  for  T.. 


Table  8.  Current  Loop  Torque 


h  (km) 

250 

500 

1000 

36000 

B  (gauss) 

4. 5  x  10'1 

4.0  x  10_1 

3.3  x  10'1 

10‘3 

Tj  (dyne  cm) 

9.0  x  106 

8.0  x  106 

6.  6  x  106 

2  x  10 

l 
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4.3  Gravitational  Torque  (Gravity  Gradient) 


Starting  with  the  structure  in  Figure  4  aligned  at  g  =  45a  to  the  local  vertical, 
we  determine  the  gravitational  or  gravity  gradient  torque  for  the  four  altitudes. 
Taking  mass  and  dimensions  from  Table  2  we  have,  since  I,x'>I,z  and  Lz'L  : 

(MI)  (MI)z 

-  6.  67  x  lO***  gm  cm2 
(MI)x  =  1.67  x  1014  gm  cm2. 

Substituting  in  Eq.  (20)  we  obtain  the  values  for  T^  given  in  Table  9. 


Table  9.  Gravitational  Torque 


h  (km)  250 

500 

1000 

36000 

Ty  (dyne  cm)  1.37  x  10 1  ^ 

1.23  x  1013 

9.  94  x  1012 

5.40  x  1010 

5.  CONCLUSIONS 

The  most  striking  result  from  the  evaluation  of  forces  acting  on  a  large  space 
structure  is  the  greatly  increased  importance  of  gravitational  torques.  At  all 
altitudes  considered,  this  mechanism  produces  a  torque  at  least  2  orders  of  mag¬ 
nitude  greater  than  any  other.  A  gravity  gradient  stabilization  mode  is,  therefore, 
unavoidable,  unless  eif  ier  a  continuously  operating  active  attitude  control  system 

is  used  or  the  structure  is  designed  to  have  nearly  equal  moments  of  inertia  about 
2 

all  three  axes. 

9  11 

In  the  first  case,  the  mass  of  the  structure  (10  -  10  gm)  means  that  large 

amounts  of  power  would  be  required  for  such  a  system  and  that  the  interaction  be¬ 
tween  the  control  system  and  the  inevitable  flexing  of  so  large  a  structure  would 

24 

pose  totally  new  problems.  The  equalization  of  moments  of  inertia  approach 
seems  more  promising  if  gravitational  alignment  is  undesirable;  however,  if  any 
other  mechanism  is  employed  to  orient  all  or  part  of  a  structure  (solar  radiation 


24.  Ginter,  S.  ,  and  llalas,  M.  (1978)  Attitude  stabilization  of  large  flexible 

spacecraft,  AIAA  Paper  No.  78-1285,  in  Proceedings  of  AIAA  Guidance 
and  Control  Conference,  Palo  \lto.  CA,  1978,  AIAA,  New  YorkT 
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pressure  using  large,  light,  movable  "sails,  "  for  example),  the  situation  is  compli¬ 
cated  by  the  need  for  sensing  and  moving  systems  that  must  operate  continuously. 

Drag  forces  are  of  little  significance  in  practice,  except  perhaps  during  the 
assembly  phase  in  low  earth  orbit.  This  is  because  of  the  need  to  ensure  an 
orbital  lifetime  on  the  order  of  tens  of  years,  which  means  that  most  structures 
will  be  placed  in  final  orbits  above  1000  km.  Since  so  many  structures  together 
with  an  ever-increasing  number  of  communications  satellites  and  relay  platforms 
will  be  placed  in  geosynchronous  orbit,  plans  must  be  made  at  this  early  stage 
concerning  the  final  disposal  of  the  structures  at  the  end  of  their  useful  life  ("solar 
sailing"  to  much  higher  orbits,  for  example).  Additionally,  contingency  plans 
must  be  made  to  minimize  the  possibility  of  the  re-entry  of  large  substructures 
that  experience  some  malfunction  during  assembly,  since  at  space  shuttle  altitudes 
orbital  lifetimes  are  no  more  than  a  few  years. 

Great  care  must  be  exercised  in  applying  the  results  of  Sections  3  and  4  to 
other  structures  and  orbits.  In  most  examples  the  orientation  of  the  structure  has 
been  chosen  to  maximize  the  individual  drag  and  torque  forces.  It  is  possible  that 
other  combinations  of  attitude,  geometry,  and  orbit  could  combine  to  reduce  the 
very  large  differences  in  the  magnitude  of  the  various  forces. 
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Appendix 

ROTATION  EFFECTS 


Additional  mechanisms  come  into  play  when  a  body  in  the  space  plasma  has 
rotation  about  some  internal  point  in  addition  to  translational  motion.  These 
mechanisms  have  been  studied  principally  for  spherical  bodies,  1  and  are  sum¬ 
marized  here  for  completeness. 

•  Aerodynamic  Torque:  A  spinning  body  will  be  decelerated  by  drag  from 

2 

the  neutral  particles  impinging  on  its  surface.  For  a  spherical  body,  the 
magnitude  of  torque  resulting  is  negligible  above  about  600  km  altitude 
when  compared  with  the  other  mechanisms  listed  here. 

•  Surface  Charge  Torque:  The  geometry  of  any  bound  charge  on  a  conducting 
body,  induced  by  translational  motion  across  a  magnetic  field,  for  example, 
maintains  a  fixed  orientation  in  that  field.  If  the  body  spins,  the  bound 
charge  flows  through  it  to  maintain  this  fixed  orientation.  The  energy  re¬ 
quired  to  overcome  the  electrical  resistance  to  this  current  flow  is  taken 

1.  Wood,  G.  P. ,  and  Hohl,  F.  (1965)  Electric  potentials,  forces  and  torques  on 

bodies  moving  through  rarefied  plasmas,  AIAA  Paper  65-628,  AIAA, 

New  York. 

2.  Davis,  A.  H. ,  and  Harris,  I.  (1961)  Interaction  of  a  charged  satellite  with 

the  ionosphere,  in  Rarefied  Gas  Dynamics,  edited  by  L.  Talbot,  p.  691, 
Academic  Press,  New  York. 
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from  the  kinetic  energy  of  the  body's  rotation,  and,  therefore,  a  decelera¬ 
tion  of  the  spin  results.  In  the  case  of  a  sphere,  however,  this  deceleration 
is  many  orders  of  magnitude  less  than  all  other  torques  due  to  rotation. 

•  Reflected  Ion  Torque:  Ions  that  strike  the  body  are  accelerated  as  they 
approach  it  by  the  potential  of  the  body  with  respect  to  the  plasma.  If  they 
are  neutralized  and  undergo  diffuse  reflection  without  significant  thermal 
accommodation,  they  have  a  net  increase  in  energy.  The  associated  mo¬ 
mentum  loss  to  the  body  applies  a  force  to  it  at  the  point  of  impact.  If  the 
ion  flow  is  not  symmetric  about  the  center  of  mass,  because  of  the  poten¬ 
tial  distribution  on  the  body,  for  example,  then  a  decelerating  torque  is 
applied.  For  a  sphere,  this  torque  is  some  3  orders  of  magnitude  smaller 
thai  the  eddy-current  and  magnetically  induced  torques. 

•  Coulomb  Torque:  Ions  that  are  not  reflected,  but  remain  on  the  body  after 
neutralization,  also  will  cause  a  torque,  if  the  moments  of  their  impact  are 
not  symmetric  about  the  body's  center  of  mass.  The  magnitude  of  this 
torque  for  a  spherical  body  is  some  3  orders  of  magnitude  less  than  the  eddy 
current  and  induced  torques  which  effectively  determine  the  changes  to  the 
initial  spin  rate  of  an  orbiting  body. 

•  Eddy-Current  Torque:  The  rotation  of  a  body  within  a  magnetic  field  gene¬ 
rates  internal  eddy  currents  that  interact  with  the  magnetic  field  via  Lorentz 
forces  to  give  decelerating  torques.  These  are  the  largest  spin-reducing 
torques  above  about  600  km  altitude.  For  the  Echo  II  satellite,  a  41  -m- 
diam,  segmented  conducting  sphere  in  a  1400 -km  circular  orbit,  the  torque 
had  a  value  on  the  order  of  10  ^  dyne  cm. 

•  Magnetically  Induced  Torque:  The  currents  flowing  between  the  ambient 
plasma  and  a  conducting  body  give  rise  to  Lorentz  forces  as  the  body  moves 
across  a  magnetic  field.  If  the  moments  of  these  forces  are  not  equally 
distributed  about  the  center  of  mass,  an  accelerating  torque  will  result, 
tending  to  spin  up  the  body.  In  the  case  of  Echo  II,  this  torque  was  due  to 
non-uniformities  in  the  electrical  resistance  of  the  various  segments  from 
which  the  sphere  was  constructed,  and  was  estimated  to  have  a  value  on  the 
order  of  5  x  10^  dyne  cm.  Because  the  induced  torque  approximately  equal¬ 
ized  the  decelerating  eddy-current  torques,  the  spin  rate  of  Echo  II 
remained  almost  constant  for  many  months. 
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